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Abstract 

We examine self-avoiding walks in dimensions 4 to 8 using high-precision Monte-Carlo 
simulations up to length N = 16384, providing the first such results in dimensions d > 4 
on which we concentrate our analysis. We analyse the scaling behaviour of the partition 
function and the statistics of nearest-neighbour contacts, as well as the average geometric size 
of the walks, and compare our results to 1/d-expansions and to excellent rigorous bounds that 
exist. In particular, we obtain precise values for the connective constants, ^5 = 8.838544(3), 
He = 10.878094(4), ^7 = 12.902817(3), fig = 14.919257(2) and give a revised estimate of 
/i4 — 6.774043(5). All of these are by at least one order of magnitude more accurate than 
those previously given (from other approaches in d > 4 and all approaches in d = 4). Our 
results are consistent with most theoretical predictions, though in d = 5 we find clear evidence 
of anomalous A^^ ^/^-corrections for the scaling of the geometric size of the walks, which we 
understand as a non-analytic correction to scaling of the general form iV^'*"''^/^ (not present 
in pure Gaussian random walks). 
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1 Introduction 



The universal properties of linear flexible polymers in a dilute solution can be modelled by the 
lattice model of self-avoiding walks (SAW), and have been studied for this purpose Q for over 
50 years |Q. Being described by the limit n ^ of the 0(n) (j)^ field theory [||, ^ places SAW 
amongst the most fundamental models in statistical physics along with the Ising (n = 1) and 
Heisenberg (n = 3) models, while its description as a non-Markovian random walk makes it of 
intense interest to mathematicians Q. SAW are also of interest to the combinatorial mathe- 
matician as a fundamental combinatorial problem. As a critical phenomenon, in the context of 
theoretical physics, the limit of the length of the walk going to oo can be thought of as equivalent 
to approaching a critical temperature (in a generating function approach the generating variable 
acts as the Boltzmann weight in an 0(n) model). Results from the associated field theory ^ 
and subsequent confirmation by numerical methods, including careful series analysis of exact enu- 
merations 1^] and statistical analysis of high precision Monte Carlo simulations indicate 
that the upper critical dimension for SAW is du = 4. Above this dimension it is expected that 
SAW behaviour is dominated by the same behaviour as occurs in Markovian random walks (in 
a renormalisation group analysis of the 0(0) cj)^ field theory both models are controlled by the 
so-called Gaussian fixed point). We expect that while dominant exponents and ratios of scaling 
amplitudes are the same as for RW the self-avoidance constraint will affect scaling amplitudes 
and corrections to scaling. In addition, apparently asymptotic l/d-expansions give reasonable 
estimates of the connective constants. On the other hand, much is known on a mathematically 
rigorous level j^, ^, thanks to the ingenuity of the lace expansion, and so the similarity of 

pure random walk (RW) and SAW behaviour can be quantified precisely to some extent. Despite 
all this non-rigorous and rigorous information several aspects of SAW in high dimensions require 
numerical investigation. First, the connective (or growth) constants, while bounded by rigorous 
arguments and estimated by l/d-expansion values (and series analysis of fairly short exact enu- 
merations), are not known precisely from Monte Carlo simulations, and so the relative value of 
different bounds is not well understood. Second, the corrections to scaling in high dimensions have 
not been investigated, and since controversies and intriguing findings occur [^] in low dimensions 
for the SAW and related models, it is of interest to clarify these in high dimensions. In any case, 
due to SAW being such a fundamental model it clearly is of interest for us to establish as complete 
a description of the behaviour of SAW as possible. 

We have simulated self-avoiding walks on the d-dimensional hypercubic lattice in dimensions 
d = 4 to d = 8 using the Pruned-Enriched Rosenbluth Method (PERM), a clever generalisation 
of a simple kinetic growth algorithm |13, 14] using a combination of enrichment and pruning 
strategies to generate walks whose weights are largely distributed around the expected peak of 
the distribution. We have utilised an implementation similar to that described in ||T^, where 
the enrichment and pruning thresholds are dynamically changed in response to the output of 
the algorithm while maintaining a constant ratio between these thresholds. For each dimension 



1 



d = 4, . . . , 8, we have generated 10^ samples of length N = 4096 and 10^ samples of length N = 
16384. While not having completely independent samples we have crudely estimated the effect 
of the dependence and so are able to give error estimates for our values. The PERM algorithm 
is particularly appropriate for high-dimensional simulations where SAW are close to RW, which 
are simply generated by a Rosenbluth-Rosenbluth approach ||l6[. While our four-dimensional 
simulations build on the careful previous work of Grassberger et al. who simulated SAW up 
to length 4000 on the four-dimensional hypercubic lattice, our higher dimensional simulations are 
without parallel. In d = 4, our simulations do not provide any further insight |15| into the 
subtle logarithmic corrections predicted in four dimensions, but simply allow us to update the 
connective and other constants in that dimension with estimates that are an order of magnitude 
better than previously obtained. In dimensions d > 4, we compare our results to the bounds of 
the lace expansion and other approaches, to the l/d-expansion, and to series analysis of exact 
enumeration data. Our main results are contained in Tables 1, 2 and 3. Apart from the precision 
of our estimates, our other contribution is to point out evidence for anomalous sub-dominant 
corrections to scaling in five dimensions (which presumably occur in higher dimensions though so 
weakly as to be not practically measurably). 

Let the number of self-avoiding walks on the lattice of interest be ctv, that is cn = \^n\ where 
f^iv is the set of all self-avoiding walks (p of length steps {N + 1 sites) with one end at some fixed 
origin. Let pat be the number of self-avoiding polygons (closed walks) of length A^. In this paper 
we consider the d-dimensional hypercubic lattice for d = 4, 5, 6, 7 and 8. We define a reduced 
free energy or rather entropy per step kat as 

nN = -^\ogCN. (1.1) 

Let {Q)n denote the simple average of any quantity Q over the ensemble set of allowed paths VLm 
of length N . Let M{(p) be the number of non-consecutive nearest-neighbour contacts (pairs of 
lattice sites occupied by the walk) for a given walk (p. We define a normalised mean number of 
contacts per step by 

"^iv = ^, (1.2) 
and a normalised fiuctuation in the number of contacts per step by 

fN = ^ . (1.3) 

For a SAW model where each configuration is weighted by a Boltzmann weig ht (say, a;*^^'^)) to 
the number of nearest-neighbour contacts (this model is known as interacting SAW or ISAW) the 
quantities niN and /at are proportional to the internal energy and specific heat in the limit lij — > 1. 
The thermodynamic limit for SAW is given by the limit A^ ^ oo so that the thermodynamic limit 
entropy per step is given by 

Koo = hm ■ (1.4) 
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Given the thermodynamic hmit exists this quantity determines the partition function asymptotics, 
i.e. ctv grows to leading order exponentially as fi'^ with the connective constant ;U = e'^°°. 

In our simulations we also calculated two measures of the walk's average size. Firstly, speci- 
fying a walk by the sequence of position vectors tq, ri, r^r the average mean-square end-to-end 
distance is 

{RDn = {{tn - ro) • (rjv - ro)) . (1.5) 

We shall use the symbol ^ to be equivalent to 

RIn = {RDn. (1.6) 

The mean-square distance of the sites occupied by the walk from the endpoint fq of the walk is 
given by 

N 

N + ~ 



i=0 

Again we define 



1 ^ 

{RDn = j—r^ Y.(^r, - ro) • (r, - fq)) . (1.7) 



RIi,N = {Rm)N- (1-8) 

Our main new results concern d > 4, so we shall now discuss the theoretical predictions for 
those dimensions. It has been proved [^] that the thermodynamic limit exists, ie. /i exists. 



Furthermore it has been proved in sufficiently high dimensions y, |10[ that 

CN = A^l^ [l + 0{n-')) (1.9) 

for any e < min {[d — 4)/2, 1) and 

Rlj^ = deN {l + 0{n~')) (1.10) 

for any e < min {[d — 4)/4, 1). 

On the other hand on a non-rigorous level from the general theory of critical phenomena 
we further expect that the numbers of walks and polygons have both analytic and non-analytic 
corrections to scaling: 

«~^,"iv-(i+^;i+^) (1.11) 

and 

with 7 — 1 = and a — 2 = —d/2 for c? > 5. From this we deduce that the entropy, mean number 
of contacts and their fluctuations scale as 

~ ^oo + ^ + ^ + J^^, (1.13) 



(1) (2) 

~ + — + ^ + (1.14) 

and 

~ /<x> + ^ + ^ + (1.15) 

respectively. The non-analytic correction to scaling exponent Ag is associated with the strongest 
non-analytic correction. One would expect further non-analytic (other A exponent-like terms) and 
analytic corrections (e.g. a N~'^ term) — see [18| for a more in-depth discussion of possible scaling 
forms in general dimensions. Our numerical studies can discern only the strongest corrections to 
scaling. Again for the geometric size of the walk one may hypothesise that 

<^~4iV-fl + | + ^) (1.16) 



with 2v = \ for d > 5. 

It may be tempting in a non-rigorous treatment of SAW above the upper critical dimension 
to implicitly assume that the non-analytic corrections to scaling either do not occur or only occur 
with the same exponents as occur in RW. However, the field theoretic description of critical phe- 
nomena above the upper critical dimension is subtle (partially due to the presence of dangerously 
irrelevant variables) and mean-field theory is unlikely to be the whole story. For example, it is 
often assumed that hyperscaling relations break down above the upper critical dimension. On the 
other hand, it is widely accepted that the relation 2 — a = dv holds for self-avoiding polygons in 
all dimensions, where 2 — a is the entropic exponent associated pN and v is the size exponent, 
in seeming contradiction. It is certainly true that dominant exponents are usually not controlled 
by the fluctuation dominated critical behaviour that give rise to hyperscaling {2 — a ^ dv for 
the Ising model for d > A) but rather mean-field energy vs entropy physics. However, it may be 
that remnants of the fluctuation driven critical behaviour still occur in high dimensions albeit 
now contributing to the corrections to scaling. In this picture the upper critical dimension du is 
the dimension below which fluctuation driven critical phenomena (characterised by hyperscaling 
relations) are dominant, while above du they are sub-dominant to mean-field criticality (fixed 
exponents). One hyperscaling relation expected to break down for d > 5 in SAW is 2A4 — ^ = dv 
with the "gap" exponent A4 associated with the "intersection" probability (see |T8| , ^ for exam- 
ple). One may be tempted to hypothesise a correction to the scaling to cat of a term jj'^ jsf^^i-d.'^-'^ _ 
Since it is accepted that A4 = 3/2 (and proved that v = 1/2) we can hypothesise a correction to 
scaling exponent arising from such a term as 

A.^^. (1.18) 

That is, in d = 5 we expect that Ag = 1/2, so corrections of order N~^/'^ as well as analytic 
corrections of order N^'^, N^"^ etc may occur. Unfortunately we are only practically able to detect 
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corrections of order N^'^ in our simulations and we have been unable to detect even the N"^/"^ in 
d = 5. We note in passing that there are predicted logarithmic corrections in d = 4 which we can 
see, to a similar extent as in Q, in our estimations of /u — in fact we have utilised this expected 
behaviour to give our refined estimation. We now comment that the non-analytic correction to 
scaling exponent hypothesised above is the same as the "crossover" exponent 0e = (f^~4)/2 of the 
Edwards model. This leads us to conjecture that while there is strictly no crossover from Gaussian 
to non-Gaussian behaviour, the excluded volume could still make itself apparent in scaling through 
a scaling function in the variable hN'^'^ , where h is the bare measure of the excluded volume. We 
concede that a further assumption about the expansion of the scaling function is needed here. In 
any case following this line of argument it is then likely that such a correction to scaling term 
occurs in other quantities such as in the scaling of the size measures. 

Assuming that the Edwards model crossover exponent provides the dominant corrections to 
scaling exponent for the size measures, -Rgjy ^'^'^ -^miV' ^^^^ gives us 

A, = = (1.19) 

Hence in d = 5 we expect that A,. = 1/2, and whenever the value of A^ coincides with an analytic 
correction to scaling (e.g. d = 6) there may also be confluent logarithms appearing. We have been 
able to successfully test for A^ = 1/2 in d = 5 (see below), and we even have some evidence of 
confluent logarithms present for d = 6. One can also predict that for large the quotients 

Bn = RI,n/RIn -B^(l + ^ + ^) (1-20) 

approach the random walk value, B^o = dm/dg = ^ with the same type of corrections as the size 
measures approach their limits in d > 5. In conclusion, from the scaling forms above we predict 
that in d = 5 we generically expect to see, within the quality of data obtained, a correction of the 
order A^~^/^ in the size measure quantities. In dimensions 7 and 8 we expect to see only the A^~^ 
corrections, while in d = 6 we may expect to see some confluent logarithmic correction term such 
as log(iV)/iV. 

Our simulations allow us to estimate //, moo, foo, dg (we certainly confirm that de = 2dm in 
all dimensions studied) and study the corrections to scaling in kn, rn^, ReN^ RmN^ Bn- Let 
us first discuss the constants as they provide the most important information contained in this 
paper. Our best estimates and various comparisons are provided in Tables ||, |^ and ^ We are 
able to compare our results to estimates and bounds from various sources. Fisher and Gaunt 
p9| ] used exact enumerations on general d-dimensional hypercubic lattices to give an asymptotic 
1/d-expansion for up to fifth order in the variable s = {2d — 1). Nemirovsky and Freed 



extended this to include de, while Ishinabe et al. |21] extended it to include moo and /oo- In all 
these expansions the error is uncontrolled, and since they are considered asymptotic expansions 
the optimal number of terms to be used to give an accurate answer varies with dimension. All 
terms have been used when applied to d > 5 since Fisher and Gaunt [19| proposed that d terms 



plus half the next should be used in general. The l/d-expansions for /i, mgo, /oo, d^, and d^ are 

^ = s-l/s-l/s"^ -11/s^ -62/s'^ + ... (1.21) 
moo = l/s + l/s2 + 7/s^ + 35/s^ + 250/s^ + ... (1.22) 
/oo = l/s + 4/s2 + 29/s^ + 152/s*^ + 752/s^ + ... (1.23) 
de = 1 + 2/s + 28/^2 + ISO/s^ + 1382/s^ + . . . (1.24) 

where s = {2d — 1) and the expansion for d^ is trivially given by half of the expansion for de- 
The specific values for d = 4, . . . ,8 are given in Tables 1, 2 and 3. On the other hand, there has 
been much effort expended to obtain rigorous and (semi-)rigorous upper and lower bounds for the 
connective constants 11, H, H, H, m in all dimensions. The lace expansion provides excellent 



lower bounds |11, ^ not only in high dimensions such as 5, 6, 7 and 8 but also in dimensions 
3 and 4: lower bounds as quoted from Hara and Slade ||ll| for dimensions 4, 5 and 6, and for 
dimensions 7 and 8 have been computed via equation (2.34) of |11], are given in Table 1. The 



best current upper bounds |J] are also listed in Tables 1, 2 and 3 — the upper bounds for 



dimensions 7 and 8 have been computed using the Maple code available at |26]. We also include 
in our tables the previous most precise estimates. In dimensions 5, 6, 7 and 8 these have been 
from series analysis of exact enumeration data |6|, |2^, 28 1, while in dimension 4 the previous best 



estimate of ^ was obtained by Grassberger P]. The previous estimates for tti-oo are from exact 



enumerations and were given by Douglas and Ishinabe |28]. 

We now turn to our evidence for the types of corrections to scaling. For the entropy, mean 
number of contacts, and fluctuations in the mean number of contacts, we find that extrapolations 
assuming dominant 1/A^ corrections produce consistent extrapolates in all dimensions d > 5 and 
all lengths > 128. Similarly for the size measure data, ]\f, n ^^'^ in dimensions d = 7 
and d = 8 the assumption of 1/A^ corrections produces consistent extrapolates for all lengths 
N > 128. Thus we conclude that 1/A^ corrections dominate in those dimensions as predicted 
above. 

The evidence for anomalous scaling is summarised in Figures |l] and for d = 5, we clearly 
detect A^^ ^/^-corrections (see Figure |l|) in the size measure data, and for d = 6 our results are 
suggestive of N^^ log A^-corrections, which produce a slightly better fit than A^^^-corrections (see 
Figure ^). We note in passing that the absence of A^^ ^/^-corrections in past extrapolations from 
exact enumeration data (see Table ^) may have affected previous estimates for A(,. 

In summary, we have presented a comprehensive study of scaling of self-avoiding walks at 
and above the upper critical dimension, testing various scaling predictions and providing precise 
estimates of associated scaling amplitudes. 
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« 



dimension 



6 



7 



estimate for ^ 



6.774043(5) 8.838544(3) 10.878094(4) 12.902817(3) 14.919257(2) 



previous estimates 
lower bound|ll] 



upper bound p5 



1 / d-expansion [ 19 1 



6.77404(4) 8.8386(8) 10.8788(9) 12.900 



6.7429 
6.8040 
6.7714 



.8285 
.8602 
.8397 



10.8740 
10.8886 
10.8800 



12.8811 
12.9081 
12.9040 



14.920 

14.9030 
14.9221 
14.9200 



Table 1: Numerical values for self-avoiding walk connective constants fj, in dimension 4 to 



Q 



dimension 


4 


5 


6 


7 


8 


estimate for moo 


0.17088(5) 


0.134576(6) 


0.106902(4) 


0.087715(2) 


0.074222(2) 


previous estimate 
1 / d-expansion 28 ] 


28] 


0.1740(15) 
0.213125 


0.141(1) 
0.142627 


0.111(1) 
0.108376 


0.0892(8) 
0.087925 


0.0744(6) 
0.074206 


estimate for /oo 


0.330(2) 


0.2324(3) 


0.1640(2) 


0.12331(7) 


0.09818(5) 


1 / d-expansion 28 ] 


0.417 


0.2362 


0.1608 


0.12114 


0.09703 



Table 2: Numerical values for normalised mean moo fluctuation f^o of nearest-neighbour 
contacts in dimension 4 to 8. 



in 



dimension 


4 


5 


6 


7 


8 


estimate for de 




1.4767(13) 


1.2940(6) 


1.2187(3) 


1.1760(2) 


previous estimate [pO|] 




1.434 


1.296 


1.222 


1.178 


1 / d-expansion ||2^ 




1.385 


1.273 


1.212 


1.174 


estimate for dm 




0.7385(6) 


0.6470(2) 


0.6094(1) 


0.5880(1) 


estimate for Boo = dm/de 


0.504(7) 


0.5001(6) 


0.5000(2) 


0.5000(1) 


0.5000(1) 



Table 3: Numerical values for distance amplitudes dm and de and their quotient B^o = dm/de in 
dimension 4 to 8. 
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0.501 



0.500 



0.499 



0.498 - 



0.497 - 



0.496 
0.10 0.00 



II- 



0.05 
7V-V2 



0.10 



Figure 
the N- 



1: Rl^/N and Br, = Rl^^^/Rl jY versus N for d = h, showing clearly the presence of 
-correction to scaling. Prom the right-hand- side plot we extrapolate Boo = 0.5001(6). 
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0.501 



0.501 



0.500 



1- 



0.500 - 



II 



■III 



0.499 



0.499 



0.498 



0.000 



0.005 



0.498 



0.010 



0.00 



0.02 
N-^logN 



0.04 



Fi gure 2: Bjsr = j^/ -^e at versus N ^ and versus N ^ log N for d = 6, showing the possible pres- 
ence of a confluent logarithm for the -correction to scaling. We extrapolate B^o = 0.5000(2). 
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